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ABSTRACT 

The field of exoplanetary science has diversified rapidly over recent years as the field has progressed 
from exoplanet detection to exoplanet characterization. For those planets known to transit, the pri- 
mary transit and secondary eclipse observations have a high yield of information regarding planetary 
structure and atmospheres. The current restriction of these information sources to short-period plan- 
ets may be abated in part through refinement of orbital parameters. This allows precision targeting 
of transit windows and phase variations which constrain the dynamics of the orbit and the geomet- 
ric albedo of the atmosphere. Here we describe the expected phase function variations at optical 
wavelengths for long-period planets, particularly those in the high-eccentricity regime and multiple 
systems in resonant and non-coplanar orbits. We apply this to the known exoplanets and discuss 
detection prospects and how observations of these signatures may be optimized by refining the orbital 
parameters. 

Subject headings: planetary systems - techniques: photometric - techniques: radial velocities 



1. INTRODUCTION 

The currently known diversity of exoplanets is greatly 
attributable to the revolution of the transit detection 
method over the past 10 years. The measurement of 
radius and hence density were the first steps from the 
results of this technique, but soon to follow were atmo- 
spheric studies from both primary transit and secondary 
eclipse. However, the unknown inclination of the plan- 
etary orbits makes this technique only applicable to a 
relatively small fraction of the known exoplanets. For 
the non-transiting planets, reflected light and phase vari- 
ations present an additional avenu e through which to 
inves t igate planetary a tmospheres (jCharbonneau et al.l 
119991 iLeigh et al.l [20031. The net result of this new in- 
formation has lead to an unprecedented ability to char- 
acterize exoplanets. 

Phase functions in the Infra-Red (IR) primarily mea- 
sure the thermal properties of the planet whereas op- 
tical measurements probe the planetary albedo. The 
relation between giant planet atmosph eres and phase 
curves h ave been descri b ed in detail by ISudarskv et all 
(|2005f) . Ilro k, Demingl (|2010h further investigate the 
time- variation of the atmospheres in the IR for eccentric 
planets using radiative transfer models. The phase vari- 
ation of our own sola r system has been investigated by 
iDvudina et al.l (|2005l ) , and it has been shown that phase 
functions can be us ed to produce longitu dinal thermal 
maps of exoplanets (jCowan fc Agoll [2008V In addition, 
ph ase curves of ter restrial planets have been considered 
bv lMallamal ([2009h . 

Phase variations of exoplanets in the IR and op- 
tical regimes have had success due to increased ac- 
cess to improved instrumentation and space-based ob- 
servatories. This has primarily been investigated for 
transiting planets since the edge-on orbital plane pro- 
duces the highest phase amplitude. Examples of ob- 
served phase v ariations in the IR (us ing Spitzer) include 
HD 189733b (jKnutson et al.1 l2009ah and HD 149026b 
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(jKnutson et alJl2009bf ). Examples in the optical inc lude 
Kepler observations of HAT-P-7b (|Welsh et al.ll2010l ) and 
phase variations det ected in the light curve of CoRoT-lb 
([Snellen et aT1l2009ft 

iKane fc von Braunl (|2008f ) and IKane fc von Braunl 
(2009) showed that planets in eccentric orbits have 

inflated trans it probabilitie s, as demonstrated by 

HP 17156b (iBarbieri et al.l [2007h and HD 80606b 
(jLaughlin et al.ll2009|) . These types of planets will pro- 
duce relatively high phase amplitudes during a brief pe- 
riod (periaston passage) of the orbit. However, phase 
variations of non-transiting planets ha ve been restricted 
to ho t Jupiters, including v And b ([Harrington et al.l 
[20061) and HD 179949b ([Cowan et al.l 120071) There 
have been searches for phase va r iation s of, for exam- 
le, HD 75289Ab dRodler et al l l2008h and r Boo b 



Charbonneau et al. 1999 URodler et al. 



20101) but no sig- 



natures were detected in either case. Exploring the at- 
mospheric properties of longer-period planets requires 
taking advantage of highly-eccentric non-transiting sys- 
tems. Thus we also constrain the inclination and hence 
the mass of such planets. 

Here we investigate the expected photometric phase 
amplitude of long-period eccentric planets and show how 
planetary orbits in resonance can result in ambiguous 
phase variation detections. We further apply this anal- 
ysis by calculating maximum flux ratios for the known 
exoplanets and considering several interesting case stud- 
ies. We also determine the effective orbital phase regime 
over which detectability is maximised and show how re- 
finement of orbital parameters can allow efficient tar- 
geted observations at these times. This study is intended 
from an observers point of view in so far as these observ- 
able variations can be reconnected back to the theoretical 
models of exoplanetary atmospheres. 

2. EXOPLANET PHASE VARIATIONS 

In this section, we establish the theoretical framework 
which will be applied in the re mainder of the paper, 
similar to the formalism used by ICollier Cameron et al] 
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Fig. 1. — Orbit of eccentric planet, showing orbit phase angles 

ure 

Q] shows a top-down view of an elliptical planetary orbit. 
The phase angle a is described by 

cos a = sin(a; + /) (1) 

where u) is the argument of periastron and / is the true 
anomaly. The phase angle is defined to be a = 0° when 
the planet is at superior conjunction ("full" phase). In 
terms of orbital parameters, this location in the orbit 
occurs when w + f = 270°. 

The flux at wavelength A incident upon the planet is 
described by 

where L+ is the luminosity of the star and r is the star- 
planet separation. This separation is given by 

r =°"- e2 ' (3) 
1 + ecos/ w 

where a is the semi-major axis and e is the orbital ec- 
centricity The geometric albedo of a planet is defined at 
a = 0° as follows 

Fr(0,A) 



(4) 



where F r is the reflected light from the planet. The plan- 
etary flux received at Earth is then 



f p (a,\) = A g (\)g(a,\)Fi(\)-Z 



(5) 



where R p is the planetary radius, d is the distance to the 
star, and g(a, A) is the phase function. Since the stellar 
flux received at Earth is 

then the flux ratio of the planet to the host star is defined 
as 

e(a ,A)EEM^) = Ag (X)g(a,X)^ (7) 

and thus contains three major components; the geomet- 
ric albedo, the phase function, and the inverse-square 
relation to the star-planet separation. Note that for a 
circular orbit, only the phase function is time dependent. 



2.1. Wavelength Dependence 

As noted in the previous section, the observed flux 
ratio from an exoplanet is wavelength dependent. In 
particular, the atmospheric composition drives the scat- 
tering properties and thus the forms of the geometric 
albedo and phase function. This dependence has b een 
considered in great detail by iSudarskv et al.l (|2005l) in 
which they construct empirical models for exoplanet at- 
mospheres and intergrate over the surface with assumed 
opacities depending upon atmospheric composition. This 
thorough analysis is not reproduced here, but we do 
use the results of their analysis for a restricted wave- 
length, particularly with regards to the albedo function 
discuessed in the following sections. Here we confine our 
study to optical wavelengths centered on 550 nm. This 
broad ly encompasses the results ob t ained by such stud - 
ies atlCollier Cameron et ail (|2002f ) . iLeigh et al.l (|2003ft . 
and lRodler et al.l (|2010t ). This also places the study near 
the peak response of the Kepler CCD, the relevance of 
which will be discussed in later sections. 

2.2. Geometric Albedo 

It has been shown through atmospheric models that 
there is a dependence of the geometric albedo of 
giant planets on th e semi-major axis of t he orbit 
(|Sudarskv et"aDl2000l [20051: iCahov et al.ll2010h We di - 
rect the reader to Figure 9 of ISudarskv et all (2005) 
which details the wavelength and star-planet separation 
dependence of the geometric albedo. Jupiter is known 
to have a visual geometric albedo of ~ 0.52. However, 
the strong irradition of the atmospheres of giant planets 
in short-period orbits results in the removal of reflective 
condensates from the upper atmospheres and thus a sig- 
nificant lowering of the geometric albedo. Observations 
of HD 209458b using the Microvari ability and Oscilla - 
tions of STars (MOST) satellite by iRowe et al.l (|2008f ) 
failed to detect phase variations and thus they were able 
to place an upper limit of A g < 0.08, s ubsequently in- 
vestig ated using model atmospheres by iBurrows et al.l 
(2008). More recent observ ations of HAT-P-7b using Ke- 
pler by IWelsh et al.l (|2010t ) revealed phase variations in 
the light curve from which they were able to deduce a 
geomtric albedo of 0.18. 

Through the examples mentioned abo ve, and the con- 
sidera tion of the theoretical models of ISudarskv et al.1 
l|2005h . we have constructed a model which approximates 
the albedo of giant planets as a function of star-planet 
separation (see Equation |3]) . This is a hyperbolic tan- 
gential function of the form 



A„ 



-(r-m 



5(e r -! + e-( r - x )) 



3 

To 



(8) 



Equation[8]is plotted in Figure[2] showing the semi-major 
axes of HAT-P-7b and Jupiter for reference. This func- 
tion well represents the rapid ri se in optical alb e do be - 
tween 0.2 and 1 AU described by ISudarskv et al.l (|2005l) . 
as well as the continued rise beyond 2 AU whereby water 
clouds begin to be present. Though broadly robust to 
encompass both the theoretical calculations and the lim- 
ited number of measured examples mentioned above, this 
empirical function does not account for planetary grav- 
ity whose variation may affect the atmospheric properties 
and thus the albedo properties. 
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star-planet separation (AU) 

Fig. 2. — Approximation of the geometric albedo distribution for 
giant planets, where A ~ 550 nm. 

2.3. Phase Function 

A planetary phase function can be considered to be 
defined by the continuous presence of an imaginary line 
conecting the center of the star and planet which is nor- 
mal to the day-night terminator of the planet. The phase 
function of a Lambert sphere assumes the atmosphere 
isotropically scatters over 2ir steradians and is described 

by 



/(a, A) 



sin a + (tt — a) cos a 



(9) 



and is thus normalized to lie between and 1. For a 
circular orbit, the phase function applied this to the flux 
ratio relation (Equation [JJ results in both a phase func- 
tion and flux ratio which are maximum at a phase angle 
of zero. Generalizing the phase angle (see Equation [1]) 
and thus the phase function for an eccentric orbit requires 
first solving Kepler's equation 

M = E-esmE (10) 

where M is the mean anomaly and E is the eccentric 
anomaly. The true anomaly is then related to the eccen- 
tric anomaly by 



cos / = 



cosi? 



1 — e cos E 



(11) 



where the true anomaly establishes the time-dependent 
variation of the phase function. 

For th e analysis performed here, w e a dopt the ap- 
proach of lCollier Cameron et al.l (|2002h and lRodler et al.1 
(|2010f ) which utilizes the empirically derived phase func- 
tion of lHiltonl (11992ft . This is based upon observations of 
Jupiter and Venus and incorporates substantially more 
back-scattering due to cloud-covering. This approach 
contains a correction to the planetary visual magnitude 
of the form 



Am{a) = 0.09(a/100°) + 2.39(a/100°) 5 
which leads to a phase function given by 
g(a) = i(r - 4Am ( Q ) 



0.65(a/100°) 3 
(12) 

(13) 



where the wavelength dependence has been removed (see 
Section l2"Tj) . This Hilton phase function is used through- 
out the remainder of this paper. 



Shown in Figure [3] are phase functions and normalized 
flux ratios for various eccentricities and orbital orienta- 
tions. Note that the maximum flux ratio does not neces- 
sarily occur at zero phase angle for a non-circular orbit. 
This is because the orbital distance is changing and in- 
deed we shall show in later sections that the star-planet 
separation component of Equation [7] becomes dominant 
for highly eccentric orbits. This time-lag between max- 
imum flux ratio and m aximum phase was also noted by 
iSudarskv et all (|2005l) . 

2.4. Orbital Inclination 

For interacting systems, many of the system param- 
eters depend on the orbi tal inclination angle , i, of the 
system (for example, see iGelino et al.l (12006)). For ex- 
oplanetary systems, given an assumed albedo, the true 
amplitude of the phase variation can be used to estimate 
the inclination angle, and therefore constrain the mass 
of the planet derived from radial velocity data. To add 
the effect of inclination angle to the phase function, the 
phase angle (Equation [lj is modified as follows: 



sin(w + /) sini 



(14) 



At 1st and 3rd quarter (a = 90° and a = 270°), the 
flux ratio is completely independent of inclination angle 
(see Figure Q}. The effect of inclination on the shape of 
the phase function is quite small and has negligible effect 
on the location of the minimum and ma ximum values of 
the flu x ratio, as shown by Figure 20 of ISudarskv et al.l 
(2005). One complicating factor in this simple inclina- 
tion consideration is that any additional light sources (i.e. 
planets) in a given system will dilute the signature from 
the dominant light-reflecting planet. The dilution of the 
signature will remain constant if the additional planets 
are in face-on (i = 0°) orbits. This will be discussed 
further in Section [3TTI 



3. MULTI-PLANET SYSTEMS 

Detection of multi-planet systems is becoming more 
frequent as we are increasingly able to probe into smaller 
mass and longer period regimes of parameter space. If 
indeed single planet systems are rare, then it is highly 
likely that the phase curve from a particular planet will 
be "contaminated" by the reflected light of other planets 
in the system. For planets which are similar in size, the 
effect of an outer planet to the combined phase curve will 
be small since (from Equation [Jj e(a, A) cx r~ 2 . Here 
we discuss the specific cases of orbital coplanarity and 
resonant orbits. 

3.1. Coplanarity 

Additional planets in a system serve to dilute the sig- 
nature from the dominant light-reflecting planet. De- 
pending upon planet formation scenarios, it cannot be 
taken for granted that planets within a system will lie in 
coplanar orbits. As shown in Section [2~4l the phase func- 
tion may be reduced in amplitude significantly for orbits 
inclined relative to the line-of-sight, to the extreme of 
eliminating a time- variable photometric signature of the 
orbit if it is face-on (i = 0°). 

Consider the case of the planetary syste m orbiting the 
star v And, which was first discovered by iButler et al.l 
(1999). A search for reflected light from the inner-most 




Fig. 3. — The phase functions (dashed line) and normalized flux ratios (solid line) for various eccentricities and periastron arguments; 
e = 0.3 and w = 0° (top-left), e = 0.3 and w = 90° (top-right), e = 0.6 and u = 0° (bottom-left), e = 0.6 and u = 90° (bottom-right). 



planet was carried out by lCollier Cameron et al.l {2002), 
but only a marginal detection was produced leaving an 
ambiguity of the result concerning the degeneracy be- 
tween the planetary radii and assumed albedos. The 
v And system is one of the few systems which has been 
monitored astrometrically as well as spectroscopically in 
order to pro vide constraints on the orbital inclinations of 
the planets (jMcArthur et al.ll201Q[ ). The inclinations of 
the outer two planets (c and d), with semi-major axes of 
0.83 AU and 2.53 AU, were measured to be ~ 8° and 24° 
respectively. The factor of 3 increase in orbital distance 
of the outer planet results in a factor of 9 less contribu- 
tion to the total planetary reflected light from the system. 
However, the relative inclinations of the planets causes 
the phase function of the outer planet to be almost 3 
times stronger than that of the inner planet. 

3.2. Resonance 

A number of systems have now been found to con- 
tain planets in eccentric orbits with some kind of reso- 
nant behaviour. An interesting example is the HD 82943 
planetary system (see Section I7T3")) . th e 2:1 resonance o f 
which has been studied in detail by iLee et all (|2006f l. 
The periodic simulataneous periastron passage of two 
planets produces a distinct signature from a phase ampli- 
tude perspective, although those moments will not nec- 
essarily be those of maximum flux rat io. In the con- 
text of radial velocity measurements, iGiuppone et al.l 



(2009) describe how resonant orbits can affect the de- 
tectability of exoplane t s. Fu rthermore, it was shown by 
lAnglada-Escude et al.l (|2010l ) that 2:1 resonant systems 
can be mis-interpreted as single-planet eccentric orbits 
when performing a fit to the radial velocity data. The 
same is true for phase curves of multi-planet systems 
where resonant orbits can effectively hide the presence 
of the outer planet in the resulting phase curve since the 
combined phase variation will be periodic with time. In 
contrast, non-resonant planets will in general produce 
a non-periodic combined signal that will resolve as two 
separate phase functions with time. 

Shown in Figure |4] are two example systems, each with 
two Jupiter radii planets. The orbits of the first system 
are in 2:1 resonance with e = 0.0 while the planets of the 
second system are in 4:1 resonance with e = 0.5. Given 
sufficient photometric precision and observing cadence, 
it may be possible to distinguish the deviant secondary 
peak of the first system and deduce the presence of the 
outer planet. However, this will be a difficult endeavour 
since the two peaks are relatively similar in amplitude. 
The second system presents an even more difficult prob- 
lem, with the combination of high eccentricity and larger 
relative semi-major axis of the outer planet leading to 
a limited observation window, higher required cadence, 
and modest phase signature from the outer planet (seen 
close to an orbital phase of 0.95 in Figure 2]). If one is 




Fig. 4. — Normalized flux ratio (solid line) for a multi-planet system in which both planets are in 2:1 resonance with e = 0.0 (left panel) 
and 4:1 resonance with e = 0.5 (right panel). The dashed line in each case represents the phase function of the outer planet. 



unable to monitor this highest peak and also discern the 
difference in amplitude with the other three peaks then 
the presence of the outer planet will remain hidden to 
the observer. Thus the derived system architecture based 
purely upon the phase variations will be incorrect. Res- 
onant systems such as these currently comprise a small 
fraction of the total number of exoplanet systems. The 
relevance of this issue will increase as radial velocity sur- 
veys sample to longer periods and as Kepler discovers 
multi-planet syst ems, of which cand idates have already 
been announced (jSteffen et al.ll2010[ ). 

4. APPLICATION TO KNOWN EXOPLANETS 

Here we apply the results of the previous sections to the 
known exoplanets. The orbital parameters of 370 plan- 
ets were extracted using the Exoplanets Data Explorer 1 . 
The data are current as of 22nd May 2010. Since the 
flux ratio is oc i?^, the unknown planetary radii for the 
non-transiting planets injects a degre e of uncertainty into 
these calcul ations. The models of iBodenheimer et all 
(|2003h and iFortnev et ail (|2007l ) show that there is a 
clear planetary radius dependence upon stellar age as 
well as incident flux and planetary composition. How- 
ever, IFortnev et al.l ()2007l ) also showed that, for a given 
planetary composition, planetary radii should not vary 
substantially between orbital radii of 0.1-2.0 AU. Since 
most of the planets we are considering here lie beyond 
0.1 AU from their parent stars (by virtue of their eccen- 
tricity) and the mass distribution peaks at one Jupiter 
mass in this region, we fix the radius for each of the plan- 
ets in this sample at one Jupiter radius, with the caveats 
mentioned above in mind. 

As shown earlier, the orbital phase at which the max- 
imum flux ratio, e max , occurs for an eccentric orbit de- 
pends upon the orbital orientation. Since Kepler's equa- 
tion is a transcendental function, the integral of Equation 
[7] must be solved numerically in order to determine the 
maximum flux ratio for each planet. These calculations 
assume an orbital inclination of i — 90° and thus the 
maximum flux ratios are upper limits in most cases, even 
though the radial velocity technique is biased towards de- 
tection of edge-on orbits since these produce larger radial 

1 http://exoplanets.org/ 



velocity semi-amplitude signatures. We have also cal- 
culated the minimum time difference in units of orbital 
phase between where maximum flux occurs and where 
the flux drops to less than 5% of the difference between 
maximum and minimum flux. This quantity is desig- 
nated At and represents the minimum time over which 
observations of maximum effectiveness can be made. 

These calculated values are plotted in FigureEJ both as 
a function of period and eccentricity. The value of At is 
~ 0.37 for all planets in circular orbits since this is where 
the phase difference for a simple cosine variation crosses 
the < 5% threshold described above. For eccentric or- 
bits At can be larger than expected, particularly where 
us ~ 270°. Even so, the distribution shown in the top- left 
panel of Figure [5] mirrors the distribution of orbital ec- 
centricities. The bottom-left panel shows that there is a 
minimum value of At that may occur for a given eccen- 
tricity, but once again we see that this can float upwards 
depending upon the value of w. The evident linear re- 
lation in log-space of the flux ratio on period shown in 
the top-right panel demonstrates that the flux ratio is 
indeed dominated by the star-planet separation as one 
would expect. However, note the significant outliers be- 
yond a period, P, of 200 days which are caused by the 
highly eccentric planets which pass through periastron 
close to a phase angle of 0° . Several of these systems are 
discussed in detail in Section [7l 

The calculated values of e max and At for ~ 70 of the 
most eccentric known exoplanets are tabulated in Table 
[TJ Of the planets represented in this table, the planet 
with the highest eccentricity, HD 80606b, is also the 
planet with the highest predicted flux ratio. This is not 
surprising considering that this planet's periastron pas- 
sage is behind the star, leading to the high secondary 
eclipse probab ility and the subsequ ent observation of 
that eclipse by Laughli n et al.l ([2009). As described ear- 
lier and demonstrated by Figure [5j the flux ratios of the 
planets in Tablc[T]arc dominated by the period and there- 
fore the semi-major axis of the orbits. 

5. ORBITAL PARAMETER REFINEMENT 

As described bv iKane et al.l <|2009f ) . the refinement of 
orbital parameters is not only an essential component 
for succussful detection of features which only appear 
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Fig. 5. — Calculated values of maximum phase amplitude e max and approximate time (in phase) between minimum and maximum 
amplitude At for the 370 known cxoplanets included in this sample. The top two panels show these calculated quantities plotted as a 
function of orbital period, and the bottom two show their variation as a function of orbital eccentricity. 



for a small fraction of the orbit, it is also achievable with 
relatively few additional radial velocity measurements. 
This is particularly true of long-period planets whose or- 
bits tend to have higher associated uncertainties and for 
which opportunities to observe at a particular place in 
the orbit are far less frequent. 

For the goal of attempting to detect a planetary tran- 
sit, it is the time of predicted transit mid-point which 
needs to be constrained. For optimal observations of 
phase variations, it is the time span during which the 
maximum change in flux ratio occurs which needs to be 
accurately determined, previously defined by the quan- 
tity At. The reason for this is because, even though the 
phase variation occurs over the entire orbit, it is assumed 
that the high cadence and precision needed will necessi- 
tate limited observing time using highly-subscribed in- 
struments (discussed further in Section [7]) . 

Using the analo gy of the transit window described by 
iKane et al.l (|2009f) , we introduce the concept of the phase 
prediction window, which is the time period during which 
a particular phase of the orbit could occur according to 
the uncertainties associated with the orbital period and 
the time of periastron passage. In Figure |6] we plot or- 
bital eccentricity as a function of the phase window for 
the 370 exoplanets for which the necessary uncertainties 
were available (358). It is clear that the planets with the 




10"* 10~ 3 0.01 0.1 1 10 100 1000 10* 



Phose Prediction Window (days) 

Fig. 6. — Orbital eccentricity as a function of phase prediction 
window for the known exoplanets. The dashed vertical line indi- 
cates a reasonable boundary beyond which observational attempts 
during a predicted phase location become difficult. 

highest eccentricities are the most difficult cases to pre- 
dict orbital phase locations. This is important because 
the most eccentric orbits tend to have much smaller At 
values and so it is essential that orbital refinement be 
used to reduce the uncertainty in the phase prediction 
and hence the size of the associated window. 
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TABLE 1 

emax AND At FOR ECCENTRIC EXOPLANETS. 
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Note. — At in units of orbital phase. 
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Change in Am (%) 

Fig. 7. — The percentage change in the values of At and 

^max as a 

function the percentage variation in the planetary visual magnitude 
correction (see Equation I12H for the planets in Table [T] 

6. ROBUSTNESS OF PHASE MODELS 

Although the phase function and albedo formulation 
adopted in this study to compute the expected phase 
variations are physically motivated, the sample of exo- 
planets with accurate measurements for these functions is 
relatively small. For example, the Hilton phase function 
(Equationll3jl is based upon the cloud maps of Venus and 
Jupiter which, although they produce similar phase func- 
tions, have their own unique cloud configurations with 
resulting slight variations in their respective phase func- 
tions. Likewise, real exoplanets could potentially exhibit 
a range of phase functions and albedo distributions be- 
yond those considered here. As is clear from Equation 
El the flux ratio is linearly dependent on the geomet- 
ric albedo. In other words, a 1% change in the value 
of A g translates into a 1% change in the value of e max . 
The geometric albedo has no effect on the time between 
minimum and maximum flux, At, since the shape of the 
phase variation is not altered. 

A change in the phase function however, will alter both 
the calculated values of At and e max . In order to quantify 
this effect, we varied the value of the planetary visual 
magnitude correction, Am, as described in Equation [T2l 
This was performed for all of the planets shown in Table 
[TJ from which the the mean of the percentage change 
in At and e ma x was calculated for those planets. The 
results of this simulation are shown in Figure [3 Even 
a substantial change in the phase function of 10% leads 
to a relatively minor change in the expected value of At 
and an even smaller impact on the expected flux ratio. 
Additionally, the magnitude of these changes are largest 
for eccenctric planets which are more sentitive to the 
location of the peaks in the phase curve. 

7. CASE STUDIES AND DETECTABILITY 

Here we present specific examples of predicted phase 
amplitudes and detectability for several of the known ex- 
oplanets. In evaluating whether these signatures are de- 
tectable or if instead calculating these signatures will re- 
main a theoretical exercise for the immediate future, con- 
sider the precision of the Hubble Space Telescope (HST) 
and the MOST satellite. Both of these telescopes have 
observed the V = 7.65 star HP 209458. The HST ob- 
servations by I Brown et al.l (|2001l ) achieved a precision of 
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1.1 x 1(T 4 and the MOST observations by ICroll et "all 
(2007) achieved a precision of 3.5 x 10 3 . In these cases, 
the necessarily high cadence resulting from the brightness 
of the host stars could be used to advantage by binning 
the data to produce higher pre cision. In addition, the el- 
lipsoidal variations detected by IWelsh et al.l (|2010h using 
Kepler data is of amplitude 3.7 x 10~ 5 . 

7.1. The HD 37605 System 

Our current knowledge of the HD 37605 system con- 
sists of a single planet in a ~ 54 day, highly eccentric 
orbit (jCochran et al.ll2004h . This places this planet close 
to the top of the list shown in Table [TJ The peak flux 
ratio from this planet is expected to be 0.56 x 10 -5 , more 
than a factor of 10 smaller than that for HAT-P-7b. Even 
so, the peak flux ratio is helped substantially by the peri- 
astron argument of u = 211° which places the periastron 
passage close to the observer-star line-of-sight on the far 
side of the star. 

The minimum time between minimum and maximum 
flux ratio is 0.04 of the orbital phase, or ~ 2.0 days. 
Conversely, the phase prediction window for this planet 
at the time of discovery (2004) was ~ 1.4 days. At the 
time of writing, this window has since grown to ~ 19 
days due to the ~ 40 orbits which have occurred since 
then. Thus, the precision of the orbital parameters re- 
quire further improvement before a robust attempt to 
only observe the At section of the orbital phase is made. 

7.2. The HAT-P-13 System 

The HAT-P-13 system ()Bakos et al.ll2009D is particu- 
larly interesting because it consists of an inner transiting 
planet (P = 2.92 days) in a circular orbit with an outer 
companion in a long-period (P = 4 28.5 days) eccentric 
orbit. Additional data acquired by I Winn et all (|2010T ) 
provide evidence for a third body in the system. The 
coplanarity of the system remains in question however 
since the refined orbital parameters for the outer planet 
have not yet resulted in a transit detection. Could this 
be resolved with observations during the predicted max- 
imum flux ratio occurance of the outer planet? This 
will be difficult for the following reasons. Firstly, the 
maximum flux ratio of the inner planet is predicted to 
be 3.11 x 10~ 5 , compared to 0.017 x 10~ 5 for the outer 
planet. Detecting the signature of the outer planet would 
be a challenging detection task even if the signature of 
the inner planet were not present. Secondly, since the in- 
clination of the outer planet is not known, the predicted 
flux ratio becomes an upper limit, potentially making 
the detection criteria even more dire. The substantiative 
dominance of the inner planet over the phase signature 
makes this far too challenging a task for any current in- 
struments. 

7.3. The HD 82943 System 

The HD 82943 system contains two pl anets in eccentri c 
orbits in a well-studied 2:1 resonance (|Lee et al.|[2006h . 
The inner and outer planets have periods of 219 and 441 
days respectively but have near identical maximum flux 
ratios of ~ 0.008 x 10~ 5 ; well beyond the detection limits 
of current instruments. The reason for the similarity in 
the peak flux ratios is due to the orbital elements in the 
system. The eccentricity of the orbits combined with the 
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Orbital Phase 

Fig. 8. — The predicted normalized flux ratio (solid line) for the 
HD 82943 system, the planets of which are in 2:1 resonance. The 
dashed line represents the phase function for the outer planet. 

resonance has led to the axes of the orbits being almost tt 
out of phase with each other. This leads to the complex 
structure for the total flux ratio variation of the system, 
shown in Figure |U Note that the periastron passage of 
the outer planet occurs behind the star, whereas it occurs 
in front of the star for the inner planet, yielding a relative 
increase in the maximum flux ratio for the outer planet, 
partially compensating for the larger semi-major axis. 

8. CONCLUSIONS 

Current generation space-missions are already detect- 
ing exoplanet phase variations in the optical (eg., Kepler) 
and the IR (eg., Spitzer). The steps these produce to- 
wards characterizing the atmospheres of these exoplanets 
are significant since they provide direct measurements 
of the atmospheric albedo and thermal properties. This 
has currently been primarily undertaken for short-period 
planets since many of these transit and produce phase va- 
rations on easily observable timescales. However, current 
radial velocity surveys are biased towards planets whose 
orbits are closer to edge-on, since larger semi-amplitude 
signatures are produced, and therefore biased towards 
planets with larger predicted phase amplitudes. 

We have shown here how time and position depen- 
dent functions for the geometric albedo and phase can 
be used to describe expected phase variations for long- 
period eccentric giant planets. There is a clear degener- 
acy with orbital inclination and resonace when consid- 
ering multi-planet systems and care must be taken to 
account for these possibilities. Applying these results 
to the known exoplanets shows that many long-period 
eccentric planets can have significant peak flux ratios, 
comparable to those of short-period planets. Addition- 
ally, the phase prediction windows of eccentric planets, 
during which observations will be optimally placed, will 
tend to be poorly contrained. The refinement of orbital 
parameters for the known exoplanets is clearly a key com- 
ponent for optimal observations of the eccentric planets 
during maximum phase amplitude. Improving the mea- 
sured orbits of long-period planets is already being un- 
dertaken by such projects as the Transit Ep hemeris Re- 
finem ent and Monitoring Survey (TERMS) (|Kane et al.1 
2009). However, most of the predicted flux ratios for 
the known planets push heavily against the boundaries 
of what is achievable with current ground and space- 



Phase Variations of Eccentric Planets 



9 



based instruments. A thorough search of all these planets 
will therefore likely need to await future generation tele- 
scopes, such as the European Extremely Large Telescope 
(E-ELT), the Thirty Meter Telescope (TMT), the Giant 
Magellan Telescope (GMT), and the James Webb Space 
Telescope (JWST). 

As more science results are released by the Kepler mis- 
sion, the study of photometric phase variations of long- 
period planets will become increasingly relevant. Not 
only is it an existing mission which has already detected 
phase variations in the light curve of HAT-P-7b, it is 
specifically looking for transiting planets where there is 
an a priori knowledge that the orbital inclination is fa- 
vorable towards maximum phase amplitude. Addition- 
ally, Kepler will eventually detect transiting long-period 
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(P > 100 days) planets where the bias will certainly be 
towards eccentric orb its since those have a high er prob- 
ability of transiting (jKane fc von Braunll2008f ). These 
future discoveries will be prime candidates to detect the 
phase variations described here. 
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